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We study the mechanism of the 'pearling' instability seen recently in experiments on lipid tubules 
under a local applied laser intensity. We argue that the correct boundary conditions are fixed chem- 
ical potentials, or surface tensions E, at the laser spot and the reservoir in contact with the tubule. 
We support this with a microscopic picture which includes the intensity profile of the laser beam, 
and show how this leads to a steady-state flow of lipid along the surface and gradients in the local 
lipid concentration and surface tension (or chemical potential). This leads to a natural explanation 
for front propagation and makes several predictions based on the tubule length. While most of the 
qualitative conclusions of previous studies remain the same, the 'ramped' control parameter (surface 
tension) implies several new qualitative results. We also explore some of the consequences of front 
propagation into a noisy (due to pre-existing thermal fluctuations) unstable medium. 



PACS: 47.20.Dr, 47.20.-k, 47.20.Hw, 82.65.Dp 
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1. INTRODUCTION 

A recent series of exciting experiments jlj] demon- 
strated a dynamic instability induced on tubules of single 
lipid bilayers by application of laser 'tweezers', whereby 
the cylindrical tubule of radius Rq modulates with a 
wavenumber given by q*Ro ~ 0.8. This instability has 
been attributed to an excess surface tension due to the 
gain in electrostatic energy when surfactant molecules, of 
higher dielectric constant than water, displace water in 
the electric field of the laser. 

The starting point for understanding this phenomenon 
is the Rayleigh instability ||-|| of a thin cylindrical 
thread of liquid with positive surface tension, whereby 
the thread can reduce its surface area at fixed volume 
by modulating and evolving towards a string of beads. 
Rayleigh calculated the preferred wavelength of a cylin- 
der of fluid in air in the inviscid and non-inertial (vis- 
cous) J3j limits, finding in the former case a characteristic 
non-zero wavenumber and in the latter case a preferred 
wavenumber of zero (or infinite wavelength). Later, To- 
motika Q calculated the instability for a viscous fluid 
surrounded by another viscous fluid, again in the non- 
inertial regime, finding that the change in boundary con- 
ditions restores a finite characteristic wavelength. See 
Olami and Granek g for a discussion of this point. The 
present problem, however, requires a much different de- 
tailed dynamical analysis which relates the flow of lipid 
molecules in the interface to the bulk flow in the sur- 
rounding fluid. An important physical ingredient is a 
new conserved quantity, the lipid on the surface. 

At present there are (at least) two theoretical treat- 



ments of the experiments of Ref. jjj . Bar-Ziv and Moses 
jl| and Nelson and co-workers || have proposed the pic- 
ture that the surface tension rapidly equilibrates every- 
where to an induced value So, and the instability pro- 
ceeds from this state. In contrast, Granek and Olami j|] 
have postulated that the correct treatment of the prob- 
lem is to impose a constant rate at which lipid molecules 
are drawn into the trap from the tubule. This loss of lipid 
is accommodated by stretching out small wavelength sur- 
face fluctuations and the result is again a uniform surface 
tension E . Goldstein, et al. (GNPS) demonstrated 
quantitatively how the equilibration of the tension in the 
tube stays 'ahead' of a shape change, so that a treat- 
ment with a constant (in time) surface tension is reason- 
able; and argued that the primary loss of area is in the 
shape instability itself, rather than through the removal 
of small-scale wrinkles. 

We propose a slightly different picture of the steady 
state before the onset of the instability, which follows 
from consideration of the experimental configuration. 
The tubules, as formed, arc several hundred microns long 
and are attached at either end to 'massive lipid glob- 
ules' [jlj of order 10/im in diameter. Hence, the tubules 
must be in contact with a reservoir which fixes the lipid 
chemical potential (or, equivalently, the surface tension). 
If we assume the system is equilibrated, it follows that 
the chemical potential for exchange between the tubule, 
reservoir, and solvent /lipid bath vanishes (8), and we may 
assume a reference chemical potential of zero or, equiva- 
lently zero surface tension. This coincides with the ex- 
perimental observation of visible thermal fluctuations on 
the tubules Q. 
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Now imagine applying a laser to the tubule. In the 
electric field of the chemical potential of a lipid molecule 
is lowered by an amount SeSDa, where D is the molecu- 
lar length, Se is the dielectric constant relative to water, 
a the area of the lipid, and £ the energy density deposited 
in the trap. Nelson et al. jjj calculated that this yields an 
energy gain per area of bilayer of So ~ 2 • 1CP 3 erg cm~ 2 , 
for a laser power of 50 mW. 

Hence there is a large reduction in the local chemical 
potential as the lipid suddenly finds it advantageous to 
move into the laser spot. The surface tension in the adja- 
cent portion of the tubule increases as lipids start to move 
out of the surface. Since the other end of the tubule is in 
contact with a reservoir at zero chemical potential, the 
final state (prohibiting, for the moment, surface undula- 
tions) must be a non-equilibrium steady state in which: 

1. Lipid is transported at constant velocity from the 
reservoir at zero chemical potential to the laser trap 
at a negative chemical potential. 

2. The chemical potential drops linearly along the 
tubule, with a gradient that balances the frictional 
drag of the bulk fluid in steady state. 

3. The local lipid concentration also varies linearly, 
since the two-dimensional lipid fluid membrane is 
compressible. 

This differs significantly from the treatments of Nelson 
et al. and Granek and Olami in that lipid must Row out 
of the anchoring globules and the chemical potential (or 
surface tension) never attains a non-zero constant over 
the duration of the experiment. In fact, prohibiting the 
shape instability, the boundary conditions specified by 
both Olami and Granek and Nelson et al. yield a tense 
final state as (a small amount of) area is drawn out of 
surface fluctuations, while the treatment of the anchoring 
globules as reservoirs yields the steady-state described 
below. W 

Several consequences follow from this observation. 
First, a chemical potential gradient suggests a mecha- 
nism for front propagation (j|J. The front starts at the 
laser spot where the surface tension is largest, and 'prop- 
agates' outward toward the anchoring globule simply be- 
cause the amplitude of the instability grows at different 
rates along the tube. Our results predict a speed of front 
propagation which is inversely proportional to the length 
of the tube, and is largest near the laser spot, decreas- 
ing to zero somewhere near the anchoring reservoirs; and 
a characteristic wavenumber which also decreases (much 
slower, see Fig. || below) away from the laser spot. 



The outline of this paper is as follows. In Section 2 we 
derive the linear concentration gradient in the absence of 
surface undulations. We predict a 'ramped', or spatially- 
varying control parameter, the effective surface tension, 
which is in fact the two-dimensional pressure whose gra- 
dient drives the flow of lipid against the viscous drag of 
the bulk fluid. In Section 3 we present a detailed micro- 
scopic picture of the uptake of surfactant by the trap, 
and argue that a competition between bending and com- 
pression energies modifies the effective surface tension of 
the trap. This leads to a prediction of a critical laser 
power for the onset of an instability. While this section 
may safely be omitted in reading this paper, it illumi- 
nates the nature of the instability by treating a realistic 
scenario for how the trap buckles to initiate flow. 

In Section 4 we discuss the implications of a slowly 
varying surface tension on the detailed calculation of 
Goldstein, et al. HQ- We also discuss front propagation 
within the picture of a surface tension gradient, which re- 
lates the problem to a large body of work on front propa- 
gation with 'ramped' parameters [jioj] . The issue of front 
propagation in this system is delicate [Q , and our results 
suggest at least two possibilities, which we briefly raise 
in this work and pose for further investigation. Depend- 
ing on whether noise (i.e. existing thermal fluctuations 
in the tubule) is present, we expect front propagation 
which is either (a) characteristic of that predicted by the 
so-called Marginal Stability Criteria (MSC) or (b) 

dominated by amplification of existing 'noise', which can 
lead to behavior reminiscent of front propagation for a 
steep enough ramp. We conclude in Section 5 by recall- 
ing the relevant timescales and frequencies, and summa- 
rizing our predictions and the differences from previous 
treatments. 



2. STEADY STATE 



In this section we calculate the steady-state configura- 
tion of a tubule under the action of an applied laser in- 
tensity, assuming the laser supplies a chemical potential 
— Eo at the laser spot. Note that this implies a reservoir 
in which to pack lipid molecules. In Section 3 we support 
this with a microscopic picture which leads to virtually 
the same results that we obtain in this section, with a 
prefactor of order one which depends on the laser shape. 
Note that there are several possible microscopic scenarios 
for initiating flow into the trap, and Section 3 addresses 
only one of these. 



1 Tf we wait long enough the trap will 'fill up' with surfac- 
tant and the chemical potential return to zero everywhere. 
However, in the present case of strong laser power the surface 
instability will have occurred by this time. See Section 3.3. 
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2.1. Equations of motion 

Changes in chemical potential S/j, are related to changes 
in surface tension by <5/i = — </>~ 1 <5£, where <f> is the 
lipid concentration (and hence <\>~ x = a is the area per 
lipid). Also, 



-P- 



(2.1) 



where p is the 2-dimensional pressure of the fluid of lipid 
molecules. 

The geometry of the system is taken as shown in Fig- 
ure 0, with the cylinder aligned parallel to the z-axis and 
r the radial coordinate. The boundary conditions are 

p(z = 0) = (reservoir) (2.2) 
p(z = L) = — E (laser spot), (2.3) 

where £o is the surface tension induced by the laser. 
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FIG. 1. Geometry of lipid tubule under localized tension. 

The Navier-Stokes and continuity equations for the 2D 
fluid of lipid molecules are 

d t (j> = -V-f>v) 
p s (d t +v-V)v = ry s V 2 v + (±rj s - 
-Vp + AT b -r. 



■7.)V(V 



(2.4) 
(2.5) 



Here n s and 7 S are 2D shear and bulk viscosities, p s is 
2D lipid mass density, and AT 6 • f is the viscous drag 
acting on the surface from the dissipative stress tensor 
T\g — \ri(y a up + V/3U a ) in the surrounding fluid. This 
flow is established in a vorticity diffusion time t v which is 
much smaller than other times in the problem. We ignore 
drag from outside the cylinder for the moment, since this 
flow essentially moves with the surfactant molecules and 
contributes relatively little to the boundary stress^ 2 ). 

With the above approximation, the boundary stress is 
given by the shear stress in the tube. For a uniform flow 
of lipid v = vz, the interior flow is Poiscuillc UM, 



u(r) 



2r 2 - i? 2 

v — m — z ' 



(2.6) 



where Ro is the tube radius and we use the no-slip bound- 
ary condition u(r = Ro) = v. Hence the stress acting on 
the surface is 



AT b -f = - 



2rjv 
Ro 



(2.7) 



Gradients of v in the z-direction change the flow pro- 
file from simple Poiseuille, but this has only a very small 
effect on the dynamics of establishing the steady state, 
primarily in the region of the laser spot, which we ignore 
for now. 

The final ingredient we need is the compressibility of 
the film, through the constitutive relation 



P = Po-X 5a 



(2.8) 



where po is the equilibrium pressure. 

Now we specialize to the problem at hand. We linearize 
the dynamic equations in v and 5<j> — <j) — <f>o, assume a ve- 
locity of the form v = v(z)z, and ignore the inertial term 
in the Navier-Stokes equation. Employing Eq. (2.8), we 
obtain 



d t 54> = -4>oV z v 

= f]V 2 z v - B<f>o X V 2 



where f] = ^rj s + j s . 

The boundary conditions are 

54>(z = 0) = 

5<j)(z = L) = <PoZ /B. 



2n 
Ri 



-v, 



(2.9) 
(2.10) 



(2.11) 
(2.12) 



where B — x 1 0q 1 is the two dimensional bulk modulus. 



2.2. Dynamics of equilibration 

Assigning 5<j)(z,t) — 54>{q,Lo)e l ^ qz ^^ and similarly for 
v(z,t), we obtain the following dispersion relation: 



iBq 2 



q 2 fj 



2v/Ro ' 



This yields ui ~ iRoBq 2 jr\ for q <g; q 
q 3> q* , where q 



(2.13) 

and to ~ B/fj for 
Hence, at long wave- 



lengths we have diffusive behavior governed by the fric- 
tion against the bulk fluid, while at short wavelengths the 
dynamics is dominated by the 2D viscosity. The crossover 
length is given by 1/q* ~ 0.1/mi, where we have taken 
■q ~ 10~ 2 g cm _1 s _1 , f] ~ 10 _6 gs _1 , and Ro ~ 0.5/zm. 
Hence in most cases of interest we are in the regime dom- 
inated by bulk fluid dissipation and may ignore f). 



' 2 'We may include this as, for example, the Stokes drag o n a 
cylinder, which increases the right hand side of Equation 
by of order 10% @. 
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We can now estimate (within linear response) the time 
to attain steady state after imposing the localized ten- 
sion by the laser as, roughly, the relaxation time of th e 
slowest mode given by the dispersion relation Eq. ( 2.13 ). 
Taking q = 2ir/L, we have 



2L 2 r] 



(27T) 2 i? B 



10" 5 s, 



(2.14) 



for B ~ 150 erg cm" 2 (12) and L ~ 100/im. This es- 
timate of B is a zero-temperature estimate and ignores 
small-scale thermal fluctuations which soften this modu- 
lus considerably [^3|,0. As we discuss in the conclusion 
and as shown in Reference Q], this effect can reduce B 
by up to three orders of magnitude, increasing t ss ac- 



cordingly, to of order 10 
vorticity diffusion time, 

' v — 

V 



■ s. We can compare this to the 



10" 



(2.15) 



where we take p = lg cm 3 . Since t ss > t v our assump- 
tion above of a uniform shear stress is reasonable. 



2.3. Steady State 

To fi nd t he steady state we equate the left hand side 
of Eq. (2.9 ) to zero, which yields a constant velocity v. 
From Eq (2.10) V z Scj) is also a constant and, applying the 
boundary conditions and the pressure constitutive equa- 
tion, Eq. (2.8), we find the following steady-state profile 



2nL 



E z 



(2.16) 

(2.17) 
(2.18) 



Thus we find that the steady state, excluding modu- 
lations of the cylinder, is a non-equilibrium steady state, 
where the lipid molecules run down a chemical potential 
gradient and the molecular spacing increases to reflect 
this changing local potential. An estimate above yields 
v ~ lfim a , where we use £o ~ 10 _3 erg cm -2 . The 
effective surface tension (or two dimensional pressure) is 
induced by the applied laser, and is non-zero only in the 
presence of flow. 



3. MICROSCOPIC PICTURE 

We have shown how lipid flow, which is a necessary 
condition for an effective non-zero surface tension far 
from the trap, follows from a boundary condition of fixed 
chemical potential at the trap. In this section we argue 
that this boundary condition requires an instability in 



the trap, and we present detailed calculations for a pos- 
sible scenario for the trap to initiate flow. We stress that 
there are several possible mechanisms, including buck- 
ling, ejection of micelles or bilayer structures, and growth 
of 'cancerous' membranes. This is surely not an exhaus- 
tive list. 



3.1. Basic Considerations 



We first note that a laser spot centered at z — L typi- 
cally has a Gaussian intensity profile [Q , which leads to 
an energy gain per area of lipid lA(z) — — Eq Q{z), with 



_ -(z-L) 2 /2A 2 



(3.1) 



The spot radius was estimated to be A ~ 0.15/im in the 
experiments of Bar-Ziv and Moses [|) . 

We can envision two scenarios after applying the laser: 



(a) 



Lipid can be sucked into the trap until the electro- 
static energy gain balances the cost of compressing 
the molecules in the bilayer. At this point the trap 
is full, flow stops, and the chemical potential (and 
surface tension) of the entire tube reverts back to 



(b) 



For a critical tension S* (Eq. |3.16 ) we expect the 
compressed section of the tubule to become unsta- 
ble with respect to buckling. For higher intensities 
the trap continues to fold to accommodate more 
lipid, initiating a flow along the tubule. This flow 
must be accompanied by a chemical potential (or 
surface tension) gradient, which drives the instabil- 
ity seen in the experiments. 



Our discussion suggests that the trap boundary condi- 
tion should contain the physics that, at a certain distance 
from the center of the trap, lipid is incorporated into folds 
to relieve the in-layer compression. A reasonable choice 
is 



a(z = L — A) = ao, 



(3.2) 



where A is a distance to be determined. This asserts that 
the area per head assumes its preferred equilibrium value 
at the point where the folding begins. 

In the next two subsections we derive the steady- 
state flow into the trap (prohibiting for the moment the 
'pearling' shape change). We first obtain a general rela- 
tion for the local area per lipid a(z), which depends on 
the trap boundary condition. We then deduce a crude 
criterion for the position z = L — A at wh ich the trap 
buckles. Applying the assumption of Eq. (3.2) at this 
position then yields the desired profile and steady-state 
flow. 
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3.2. Detailed Steady State 

Let us examine the steady state. The continuity equa- 
tion, Eq. (2.4), yields the condition 



v(z) — Ca{z), 



(3.3) 



where C is a constant to be determined. Hence the 
steady-state Navier-Stokes equation, ignoring the 2D vis- 
cosity, becomes 



() ^ +f iV 2fl (z)-VW(z). 



R 



(3.4) 



Distance along tubule 



5a area per lipid 
U(z) Laser potential \\ 
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FIG. 2. Profiles for steady-state deviation Sa in the area 
per lipid and the potential set up by the laser after buck- 
ling has occurred. The dotted line shows the continuation of 
Sa (Eq. 3.5) into the buckled trap region beyond the point 
2 = L — A, and should actually be replaced with buckled 
membrane at Sa — 0. The boundary conditions are: (1) 
Sa = at the reservoir, due to fixe d ch emical potential; and 
(2) Sa(z = L — A) = (see Eq. 3^). In steady state Sa 



increases linearly away from the reservoir because the mem- 
brane is under tension, which varies linearly to counteract the 
viscous drag due to lipid flow. 

The solution of this equation with the boundary condi- 
tion a(0) = ao, which follows from contact with a reser- 
voir at ambient pressure, is 



Sa 

A = 

7 



A :: 



2rjxC 

Ro 
Eo 
B ' 



(3.5) 
(3.6) 
(3.7) 



where A is determined by the boundary condition at the 
trap. Note that 7 ~ 10~ 5 , since Eo 
B ~ 150 erg cm" 2 . Using Eq. (3.2), we find 



10 3 erg cm 2 and 



Jo A ds((s) 



(3.8) 



where we have expanded for small A and made use of 
C(0) = far from the trap. For all practical purposes 
\ = -j((L-A)/L. 

Since XL <C 1 for most cases of experimental interest, 
the calculation of Section 2 applies in the region outside 
the trap (see Fig. ||) with the surface tension replaced by 



E = EoC(£-A). 



(3.9) 



For a Gaussian shape ((L - A) < 1 [£(£ - A) ~ 0.31 in 
Fig. U , so the modification to the naive boundary condi- 
tion is rather minimal. 



3.3. Trap Boundary Condition 

To complete our discussion we estimate the stability 
against folding inside the trap. This deter mine s the po- 
sition A at which the boundary condition (3.2) applies, 
as well as a critical effective surface tension parameter 
(or laser intensity) E* at which the system initiates flow. 
We imagine that the system has attained a steady state 
in the absence of buckling and flow, given by Eq. ( |3.5| ) 
with A = 0: 



a{z) = a [1 -7C0*)] ■ 



(3.10) 



Here 7 £ (z) is a measure of the compression. 

Rather than calculating the stability of a patch with 
a non-uniform area per head a(z), we calculate the sta- 
bility against buckling of a patch with uniform a = tpao, 
with if] = 1 — e, and use the resulting critical strain e* to 
determine A through 



7 C(L-A)^e* 



(3.11) 



We consider perturbations R(z) = Ro(l + u(z)) which 
preserve the volume of the fluid. This constraint yields 
the condition J[u(z) 2 + 2R u(z)} = The free en- 
ergy, which includes in-plane compression and bending, 

is 0) 



2F = I d\ 



B 



Tpa — ao 
ao 



(3.12) 



where H is the mean curvature, ao 
ipd^. For the perturbation above M, 



dzRo, and a 



' 3 'We do not include a term involving interaction with the 
laser, since we are interested in an instability at fixed particle 
number on the membrane. 
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d\ = dz(R a + u(z))y/l + u'(z) 2 

{l + u'iz) 2 )- 1 ' 2 u"{z) 



H 



Ro + u{z) 



(l + u'(z) 2 ) 3 / 2, 



(3.13) 
(3.14) 



where u'{z) = du/dz. 

To quadratic order in u(z) — 'Yli q u{q)e lqz the energy 
per unit area A becomes || 



3k „ „ 
-TT - r oBR, e 

Kq 



-f [ 3BR e + A 
V Ro 



(3.15) 



where q = qRo- The vanishing of the term in square 
brackets defines e* , the minimum strain above which this 
energy is unstable to undulations. Combining this con- 
dition with our estimate for how e{z) varies away from 
the trap, Eq. (3.11), we find the following relation which 
determines A: 



1 3 + 2q 4 



A) 



3q 2 



(3.16) 
(3.17) 



where q ~ ttRq/A. For a > er* the tubule should buckle 
inside the trap. Fig. || shows g(q) for the Gaussian laser 
spot. 
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FIG. 3. Buckling criterion inside the laser trap for a laser 
spot size A = 0.15 /wn and tubule radius Ro = 0.5 pm. 
Here, we find a* ~ 21, A ~ 0.9i? , £(L - A) ~ 0.31, and 
E ~ 0.97E . 

Our criterion depends on the trap shape, but not the 
compression modulus B. This happens because, while 
the strain induced in packing lipid in the trap varies as 
1/B, the critical strain at which buckling occurs is also 
inversely proportional to B, and the B dependence can- 
cels out. In fact, the same order of magnitude estimate 
emerges from a comparison of bending and effective sur- 
face tension. What we have gained, however, is a picture 
of the forces at play that induce the buckling. 



While our estimate apparently fails to predict buck- 
ling for typical values, relaxing a few approximations we 
have made should change this picture. First, we have as- 
sumed uniform coverage at the lowest value of the actual 
nonuniform coverage in the quiescent trap. Second, we 
have assumed an axisymmetric deformation. This is ob- 
viously not the case if the laser spot size is smaller that 
the tubule diameter. In addition, the volume constraint 
must be handled differently. Removing these approxima- 
tions should, in both cases, result in a smaller £*. For 
example, in the limit of small trap sizes we can ignore 
curvature and ask about the stability of a flat interface 
against buckling, for which the criterion above becomes 



S* = 



-2 



t(L-A) 



(3.18) 



demanding that the characteristic buckling wavevector 
q be roughly the inverse of the trap size A determines 
the critical trap size and intensity. This is essentially the 
same estimate given by Bar-Ziv, Frisch, and Moses in a 
somewhat different context lITcJ. 



4. DISPERSION RELATION AND FRONT 
PROPAGATION 

4.1. Dispersion Relation 

We have shown thus far that, under steady-state condi- 
tions before any macroscopic shape instability occurs, the 
proper boundary conditions imply a surface tension gra- 
dient along the tubule which supports lipid flow. We now 
turn to the effects of this gradient. Rather than repeating 
the the analysis of GNPS with a non-uniform surface 
tension, we note that the characteristic wavenumber at 
which the instability occurs is typically q*Ro — 0.8. Since 
Rq <C L, we suspect that the assumption of a locally con- 
stant surface tension along the tubule is a good first step. 
This allows us to transcribe the results of Ref . . 

The primary result of interest is the growth rate iv(q) 
of an undulation u(q, t), where q denotes a Fourier mode 
along the tubule. This frequency is defined through 



d 



(4.1) 



where the convective term arises because the lipids have 
an average velocity. 

In the original instability calculation presented by 
Rayleigh 0, and as has been emphasized in Refs. 
SUIJI, the structure of u(q) is as follows: 



tu(q) = $(q)T(qRo). 



(4.2) 



The function T(qR ) is determined by the energetics of 
the problem and, in our case, is non-zero for HR 2 /n 
greater than a critical value of order unity MM- In the 



G 



Rayleigh case the instability occurs for S > 0. The func- 
tion &(q) is determined by the dynamics of the problem, 
and it is here that much of the interesting and surprising 
physics lies. Energetics tells us that the most unstable 
modes are at low q, where undulations are the least 'vio- 
lent', while dynamic considerations severely penalize the 
growth of modes in the limit q — ► 0. 

Goldstein et al. Q| calculated u>(q) for a uniform sur- 
face tension S, including the effects of bending as well 
as friction between the two bilayer leaves. Changing the 
boundary conditions of their work to allow for flux from 
the reservoir adds a convective term to the dynamics, 
and, aside from our local approximation above, changes 
nothing else. A plot of w(q) is shown in Figure [| for 
several values of the surface tension a = J^Rq/k, with 
values for the bilayer friction and bulk modulus taken as 
in GNPS. 
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FIG. 4. Dispersion relation ui(k), with parameters as cal- 
culated by GNPS, j3 = 3.5,e = 0.5. Reprinted from Ref. [7]. 
Alternatively, we may consider this as a plot of ui(k, z) for 
a = 50 and z — L, L/5, L/10, L/25, and L/50. Here, k is in 
units of 2n/Ro- 



Since a is z-dependent (Eq. the growth rate w* 

and wavenumber q* of the fastest growing mode are z- 
dependent, and are greatest near the laser spot, as in 
Fig. [| A single Fourier modulation has, locally, the fol- 
lowing form: Q| 



U^Z, t) = Uq6 



iq(z — vt)-\-uj(q,z)t 



(4.3) 



where u)(q, z) is the function plotted as to(q) in Fig. 5 of 
Reference J7j and is reproduced here in Fig. |j. The z- 
dependence comes through the z-dependence of a. Note 
that we rely strongly on the condition qL 3> 1 [Note that 
in the experiments [Q with, say, L ~ 200 /zm and a diam- 
eter of 1 //m, this condition is easily satisfied, qL ~ 160.] 
Given a dispersion relation which depends on position, 
there are several immediate naive predictions: The lo- 
cal wavenumber and apparent growth rate of the pat- 



tern should decrease as the anchoring globules are ap- 
proached, with the instability vanishing at a point close 
to the reservoir where the induced surface tension is not 
strong enough to overcome the barrier due to bending, 
S* ~ k/Rq. Hence, in the experiments with a ~ 20 M, 
this occurs at I /20th of the distance from the anchoring 
globule to the laser trap. At points closer to the glob- 
ule any undulation is a decaying remnant of the pattern 
developed closer to the trap. 
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FIG. 5. Frequency lu* and dimensionless wavenumber k* 
of fastest growing mode as a function of position along the 
tubule, obtained from Fig. H by taking a{z — L) — 20. 



4.2. Front Propagation 

We now confront the issues of front propagation and 
wavelength selection. GNPS argued that the Marginal 
Stability Criterion hypothesis provides a reasonable esti- 
mate for both the propagation speed vj and the selected 
wavenumber k* JtJ. A naive extension of this calcula- 
tion, again assuming a local dispersion relation, predicts 
a spatially varying front speed and selected wavenum- 
ber. However, this prediction and the qualitative picture 
relies on two assumptions: (1) the existence of a prop- 
agating front, and (2) the absence of noise — i.e., that 
the propagation occurs into a uniform, unstable medium 
with no thermal fluctuations. The latter assumption is 
obviously not correct in detail, as thermal fluctuations, 
including modes of wavelengths comparable to the most 
unstable modes, are apparent in the experiments prior to 
the onset of the instability. Here we discuss these issues 
in the context of a spatially varying control parameter. 
Because a full treatment of the problem does not yet 
exist in the literature and is beyond the scope of this pa- 
per, we limit ourselves in this work to some numerical 
experiments and suggestions which, we hope, will stim- 
ulate further research on both this specific problem and 
the general aspects of front propagation into spatially- 
varying media in the presence of noise. For the rest of 
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this paper we refer to 'noise' as a set of random initial 
conditions which obey a Boltzmann distribution, and do 
not consider temporally fluctuating noise. 

Kramer et al. jit]] , followed by others ||l7|-|l9|l , showed 
that, in the presence of a 'ramped' control parameter 
that becomes subcritical at some point, as happens near 
the reservoir, the uniquely-selected wavenumber need not 
correspond to that determined by the MSC. This effect 
is expected to take precedence over the MSC-determined 
wavenumber at times after which non-linearities become 
important and the 'phase' of the pattern has time to dif- 
fuse of order the system size. However, we are concerned 
with the more fundamental issue of the existence of a 
propagating front. 

In the presence of 'noise' which, in the present ex- 
periment, corresponds to existing thermal fluctuations 
around the reference smooth cylindrical state, a propa- 
gating front can be expected to exist for times less than 
the characteristic growth times of existing fluctuations 
in the vicinity of the most unstable mode. Hence, given 
a quench into an unstable 'ramped' state with an ini- 
tial perturbation near the laser spot, propagation away 
from the perturbation occurs for an initial period of time, 
followed by rapid growth all along the cylinder as the ini- 
tial conditions ('noise' of unstable wavelengths) are am- 
plified to visible length scales. An initial perturbation 
near the laser spot is natural because, in practice, the 
laser spot diameter is smaller than the tubule diameter 
and a 'pinching' effect results whereby surfactant flows 
around the circumference of the tubule (as well as along 
the cylinder diameter) to fill the trap. 

The effect of a 'ramp' in the control parameter should 
be most dramatic after the noise overwhelms the front 
propagation: for a flat control parameter (no ramp) the 
noise grows randomly everywhere, and the 'front' should 
break down when the noise has grown to visible am- 
plitudes. However, for a steep enough ramp the non- 
uniform amplification of the noise could resemble front 
propagation. 

To check these conjectures we have employed a simple 
caricature of the tubule dynamics, specified by 

d 

(-g- t + vd x )u{x, t) = [a{x)2kldl + dt) u{x, t) - gu{x, tf , 

(4.4) 

where a{x) is a spatially varying control parameter cho- 
sen to mimic the dispersion relation and position depen- 
dencies in Figs. ^ and S. A choice which gives reasonable 
qualitative agreement is 



physics of front propagation into an unstable inhomo- 
geneous medium, as occurs in these experiments. For 
Fig. ||, x ~ 0.05L and a — 1/8 are reasonable. Fig. ^ 
shows the local dispersion relation lo*(x) for various a. 
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Local dispersion relation u* for trial form, 



Eq. (4.5), for a = 1/2, 1/4, 1/8. The trap is at z/L = 1. 



We have chosen the simplest possible non-linearity to 
stabilize the system. One may choose a more physical 
non-linearity such as the driving force arising from terms 
of higher than quadratic order in the mean curvature 
[Q, but our purpose here is primarily to illustrate some 
qualitative behavior of a front propagating into a noisy, 
non-uniform media. 



Figs. fj] and || show s nap shots in the evolution of the 
system, given by Eq. (4^), for an initial perturbation 
at the trap of 1% of the final amplitude (as determined 
by the non-linear term g) and an initial condition (or 
'noise') which is taken to be a superposition of 300 har- 
monics weighted with a Boltzmann weight corresponding 
to a non-zero surface tension (i.e. with an energy pro- 
portional to q 2 ). We have chosen a system size of 150 
wavelengths, and arbitrarily chosen the vertical scale to 
fill the figures. 



a(x) 



xo 



a+l 



(x - X ) 



(4.5) 



where xq is the point at which the system is absolutely 
unstable. We emphasize that this is a toy model whose 
details do not correspond to the Bar-Ziv et al. ex- 
periments, but which we believe contains the essential 



The general features are as described above: a front 
'propagates' for an initial time from the initial perturba- 
tion, after which the 'noise' takes over and a very irreg- 
ular growth quickly overtakes the system. 
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FIG. 7. Consecutive snapshots of u(x, i) as determined by 
Eq. (4.4). The conditions are: convection v = 0; initial per- 
turbation at the laser spot u(L,t = 0) = 10 _1 u(L, oo); ramp 
parameter a = 1/2; initial mean noise amplitude a fraction 
10 -6 of the final amplitude. Time intervals are every 80 time 
steps. The vertical scale is chosen to fit the amplitude of the 
undulation, and as such is a different scale than the horizontal 
scale. All snapshots have the same vertical scale, (a) and (b) 
show a propagating front; (c) and (d) show the acceleration 
due to amplification of the noise, and (e) shows a return to 
apparent propagation. 

The steeper ramp (a = 1/2, Fig. 0) has a better de- 
fined growth in the 'noise' regime, and could almost be 
called a 'front'. In contrast, the growth into the shal- 
low ramp (a = 1/8, Fig. ||) is more ragged and it would 
be charitable to call this a front. The shallower ramp 
has a very slightly faster propagation speed in the ini- 
tial regime, and an obviously faster 'propagation' speed 
in the noise-dominated regime. Both of these behaviors 
may be traced to the faster overall growth rate for a shal- 
lower ramp (where a larger fraction of the tubule is more 
unstable). 



if' 
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reservoir laser spot 

FIG. 8. Same parameters and time steps as in Fig. |^, ex- 
cept with ramp parameter a — 1/8. 

Fig. ^| shows the the results of fixing the ramp and 
varying the noise amplitude. For a noisier system the ef- 
fective propagation speed in the noise-dominated regime 
is faster, and the breakdown of the simple propagating 
regime occurs earlier. The propagation velocity before 
the noise takes over is independent of the noise ampli- 
tude. The delay before noise-dominance increases log- 
arithmically with increasing noise amplitude, consistent 
with the simple argument that the propagating solution 
exists until the noise has grown to a given amplitude, 
since this initial growth is exponential. 

To summarize, we have performed exploratory numer- 
ical calculations to investigate some of the consequences 
of a ramped control parameter, with an initial localized 
perturbation and initial global 'noise' for an initial con- 
dition, finding: 

1. At early times a front propagates away from the lo- 
calized perturbation. We find a dimensionless front 
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velocity of vu)(ko)/ko = 3.3, while the Marginal 
Stability Criteria |J predict 4.6. A similar agree- 
ment was found in the simulations of GNPS fit!. 



5.1. Physical Picture 
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FIG. 9. Front velocity vs. time, averaged over 100 realiza- 
tions of initial noise with weights described in the text. The 
velocity was measured by tracking the leading edge of the en- 
velope of wavelets, and is plotted in units of the characteristic 
velocity qj(fco)/fc o, w here w(fc) = 2 k\ k 2 — k 4 is the dispersion 



relation for Eq. (4.4) for a(x) — 1. Ramp parameter a = 1/8. 
The noise is given as a fraction of the final amplitude. 

2. In this first regime the propagation velocity is in- 
dependent of noise amplitude. 

3. After a time, which may be taken to be the time 
required for the 'noise' to grow to an observable 
amplitude, the unstable pattern rapidly develops 
everywhere on the tubule. 

4. The speed and qualitative character of this growth 
depend on the noise and ramp characteristics. 
The growth is faster for a shallower ramp and/or 
stronger noise, and looks reminiscent of a front for 
a steeper ramp; and in all cases is much faster than 
the simple front propagation from the initial local- 
ized perturbation at the laser trap. 

We have also performed calculations with no initial 
localized perturbation but, as this is probably not physi- 
cally relevant, we do not report the results here. This ini- 
tial study raises several questions which we feel are worth 
pursuing. In cases where the noise is weak enough and 
an apparent front exists, can this be understood quanti- 
tatively in terms of the gradient in the control parame- 
ter, and how does this relate to previous investigations 
of 'ramped' control parameters p0| , p^ l9[? Can a ramp 
stabilize an advancing front? 



5. CONCLUSION 



We have given the following picture of the action of 
lipid tubules upon the application of laser tweezers. In 
the absence of buckling, the laser induces a local com- 
pression of lipid molecules in the laser spot. This takes 
place in a time of order t ss ~ 10~ 5 s. A sufficiently large 
laser intensity induces a local buckling of the membrane 
in the trap, which initiates flow down the tubule from 
the reservoir. We do not have an estimate of the delay 
time for the instability in the trap. In the absence of 
undulations outside the trap, this flow would build up 
to a steady-state value v ~ 1 /mi s _1 in a time order t ss . 
The physics of this flow is a balance between drag against 
the bulk fluid and a force due to the gradient imposed 
by the chemical potential drop between the reservoir and 
the trap. 

Given the steady-state tension profile within the mem- 
brane and the reasonable assumption that the gradient 
occurs over a length (L) much larger than the critical 
wavelength (~ R ), the analysis of GNPS leads to a 
Rayleigh-like instability to undulations. This instability 
initiates near the laser spot where the chemical potential 
(or surface tension) is lowest, and propagates away from 
the spot to a point along the tubule at which the local sur- 
face tension falls below the critical tension E cr = k/Rq 
which characterizes the instability. Typical growth fre- 
quencies are uj ~ 25 s -1 0, which corresponds to times 
t w ~ 10 _1 s. We note that the experiments find a signif- 
icant time delay of order seconds before the instability, 
pf]| , a still- unexplained observation. 

The estimate above for t ss assumes a bare 2D compres- 
sion modulus B, while GNPS (see |0) pointed out that 
B undergoes significant softening at the lengthscale of 
the tubule, due to the thermal fluctuations at low surface 
tension, and estimated a decrease of up to three orders 
of magnitude. This, correspondingly, would increase the 
value of t ss to 10~ 2 s, so that an accurate quantitative 
calculation must include the dynamics of the increase in 
surface tension. This leads to, effectively, a smaller ap- 
plied tension a and hence a slower propagation speed. 
This spirit was followed in the approach of Granek and 
Olami g. 

Our interpretation assumes that the trap accommo- 
dates material by folding, or some other means. Our 
analysis suggests that the proper boundary condition 
should be a fixed surface tension S < Eo at the la ser spot, 
where the laser shape determines E through Eq. (3.9) and 
A through Eqs.( p7ll| , p6| ). 

This reservoir picture suggests that, upon turning off 
the laser, the system can revert to the original tubule by 
unfolding or, if severe topological changes have occurred 
(by, for example, budding in the laser spot or the cre- 
ation of metastable 'pearls' as seen in the experiments), 
attain some other long-lived metastable state. 
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5.2. Discussion 

In this work we have made several assumptions. The 
assumption that we can treat the anchoring globules as 
reservoirs presupposes that any damping processes re- 
tarding the transfer of lipid to and from the globules 
is negligible relative to other dynamical processes. We 
expect this to arise from the same source as the two- 
dimensional surface viscosity, which we have argued in 
Section 2.2 to be negligible. We have given a simplified 
picture of the scenario of trap buckling, where we take a 
single characteristic buckling wavevector, and treat the 
trap as uniform. This 'single-mode' approximation may 
be naive, and preliminary calculations suggest that the 
system is in fact less stable than this simple analysis 
would suggest |Q . There are also several other possible 
modes of instability which we have only mentioned but 
which could certainly play a role. We have also specified 
a boundary condition at the trap whereby the the lipid 
relaxes to its preferred area per head group, which seems 
reasonable but is not otherwise justified. Finally, we have 
made a local approximation for the variation of the sur- 
face tension so that we may use the results of GNPS. 
This applies for sufficiently long tubules, L/R 3> 1. 

Front propagation and the detailed effects of propa- 
gating into a spatially-varying medium have only been 
touched upon in our numerical treatment. This study 
still leaves much to be resolved; one important question 
is how to accurately treat the non-linear regime. This has 
been treated in different ways by Olami and Granek |^] 
who considered the non-linear effect of removing surfac- 
tant from the membrane in the absence of a gradient, and 
by Goldstein et al. JtJ, who added the correct non-linear 
terms in the bending energy to examine the propagation 
of the pearls. 

The primary new ingredients in our theory are (1) our 
treatment of the anchoring lipid globules as reservoirs 
and (2) our exploratory treatment of the role of pre- 
existing thermal fluctuations (noise) in determining the 
'front-like' characteristics of the instability. Both Nelson 
and co-workers [BQ and Olami and Granek || 'turn off' 
the reservoir. In the latter case material is drawn out of 
the existing thermal fluctuations, while Nelson and co- 
workers attribute the area change primarily to the shape 
instability itself. Olami and Granek impose a constant 
flux boundary condition at the trap, while Nelson and 
co-workers impose a fixed chemical potential — £ at the 
trap which, fairly rapidly, reduces the chemical potential 
everywhere to — £. Our picture essentially gives the same 
boundary condition at the trap, but the treatment of the 
globule as a reservoir changes the qualitative picture dra- 
matically. 

Our theory differs from previous theories in several re- 
spects, and there are many consequences which may be 
checked experimentally. Obviously, we expect flow when 
an instability develops. This could be visualized by, for 
example, fluorescence spectroscopy with a very dilute 



fraction of labelled lipids. The inhomogeneous surface 
tension implies that the local dispersion relation is also 
spatially-dependent, as in Fig @, which implies that the 
velocity of front propagation v f (which is proportional to 
w * @) an d characteristic wavenumbers should decrease 
farther away from the laser spot. Note that the character- 
istic wavelength changes very gradually compared to the 
speed of propagation, and as such would be more difficult 
to detect. It would also be interesting to see, experimen- 
tally, whether fluctuations are actually strong enough to 
destroy the front-like character, or whether two charac- 
teristic regimes exist in the experiments, as indicated in 
Fig. ^|. Finally, we mention that the opportunity of using 
laser pulses to control flow within lipid and other systems 
presents amusing possibilities and applications. 
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